*This file produces the graphs for Figure A5

preserve


gen offrezpop74 = offrespop
replace offrezpop74 = . if year !=1974

gen agpct74 = agpct
replace agpct74 = . if year !=1974

gen devpct74 = devpct
replace devpct74 = . if year !=1974

gen n = 1

collapse (mean) STARTDUM  tMean pptMean  resolution SoilStd ElevMean Ruggedness SoilMean StreamDist casino_year credit_yr agpct74 devpct74 offrezpop74 adj_start (max) post (sum) n, by(rezname)

ren post Treated
gen lnpop = ln(offrez)


*Panel A: 1974 Agriculture
twoway lfitci  agpct adj if Treated ==0, lcolor(maroon) lpattern(dash) alwidth(none) fcolor(gs13)       fintensity(inten50) || ///
lfitci agpct adj if Treated ==1, lcolor(maroon) alwidth(none) fcolor(gs13) || ///
scatter agpct adj if Treated == 0, mcolor(navy) msymbol(circle_hollow) || ///
scatter agpct adj if Treated == 1, mcolor(navy) msymbol(circle)  ///
graphr(color(white)) legend( order( 5  "Untreated" 6 "Treated" 2 "Linear Fit, Unreated" 4 "Linear Fit, Treated" ) cols(2))  ///
xtitle("Adjudication Start Date") ///
ytitle("Mean % Agriculture in 1974")


*Panel B: 1974 Development
twoway lfitci  devp adj if Treated ==0, lcolor(maroon) lpattern(dash) alwidth(none) fcolor(gs13)       fintensity(inten50) || ///
lfitci devp adj if Treated ==1, lcolor(maroon) alwidth(none) fcolor(gs13) || ///
scatter devp adj if Treated == 0, mcolor(navy) msymbol(circle_hollow) || ///
scatter devp adj if Treated == 1, mcolor(navy) msymbol(circle)  ///
graphr(color(white)) legend( order( 5  "Untreated" 6 "Treated" 2 "Linear Fit, Unreated" 4 "Linear Fit, Treated" ) cols(2))  ///
xtitle("Adjudication Start Date") ///
ytitle("Mean % Developed in 1974")


*Panel C: 1974 Population
twoway lfitci  lnpop adj if Treated ==0, lcolor(maroon) lpattern(dash) alwidth(none) fcolor(gs13)       fintensity(inten50) || ///
lfitci lnpop adj if Treated ==1, lcolor(maroon) alwidth(none) fcolor(gs13) || ///
scatter lnpop adj if Treated == 0, mcolor(navy) msymbol(circle_hollow) || ///
scatter lnpop adj if Treated == 1, mcolor(navy) msymbol(circle)  ///
graphr(color(white)) legend( order( 5  "Untreated" 6 "Treated" 2 "Linear Fit, Unreated" 4 "Linear Fit, Treated" ) cols(2))  ///
xtitle("Adjudication Start Date") ///
ytitle("ln(Off-Reservation Population in 1974)")



*Panel D: Date of First Casino
twoway lfitci  casino_year adj if Treated ==0, lcolor(maroon) lpattern(dash) alwidth(none) fcolor(gs13)       fintensity(inten50) || ///
lfitci casino_year adj if Treated ==1, lcolor(maroon) alwidth(none) fcolor(gs13) || ///
scatter casino_year adj if Treated == 0, mcolor(navy) msymbol(circle_hollow) || ///
scatter casino_year adj if Treated == 1, mcolor(navy) msymbol(circle)  ///
graphr(color(white)) legend( order( 5  "Untreated" 6 "Treated" 2 "Linear Fit, Unreated" 4 "Linear Fit, Treated" ) cols(2))  ///
xtitle("Adjudication Start Date") ///
ytitle("Date of First Casino")


*Panel E: Date of First Banking Access
twoway lfitci  credit_yr adj if Treated ==0, lcolor(maroon) lpattern(dash) alwidth(none) fcolor(gs13)       fintensity(inten50) || ///
lfitci credit_yr adj if Treated ==1, lcolor(maroon) alwidth(none) fcolor(gs13) || ///
scatter credit_yr adj if Treated == 0, mcolor(navy) msymbol(circle_hollow) || ///
scatter credit_yr adj if Treated == 1, mcolor(navy) msymbol(circle)  ///
graphr(color(white)) legend( order( 5  "Untreated" 6 "Treated" 2 "Linear Fit, Unreated" 4 "Linear Fit, Treated" ) cols(2))  ///
xtitle("Adjudication Start Date") ///
ytitle("Date of First Banking Access")


*Panel F: Stream Distance
twoway lfitci  StreamDist adj if Treated ==0, lcolor(maroon) lpattern(dash) alwidth(none) fcolor(gs13)       fintensity(inten50) || ///
lfitci StreamDist adj if Treated ==1, lcolor(maroon) alwidth(none) fcolor(gs13) || ///
scatter StreamDist adj if Treated == 0, mcolor(navy) msymbol(circle_hollow) || ///
scatter StreamDist adj if Treated == 1, mcolor(navy) msymbol(circle)  ///
graphr(color(white)) legend( order( 5  "Untreated" 6 "Treated" 2 "Linear Fit, Unreated" 4 "Linear Fit, Treated" ) cols(2))  ///
xtitle("Adjudication Start Date") ///
ytitle("Mean Distance to Nearest Stream (meters)")

 

*Panel G: Precipitation 
twoway lfitci  pptMean adj if Treated ==0, lcolor(maroon) lpattern(dash) alwidth(none) fcolor(gs13)       fintensity(inten50) || ///
lfitci pptMean adj if Treated ==1, lcolor(maroon) alwidth(none) fcolor(gs13) || ///
scatter pptMean adj if Treated == 0, mcolor(navy) msymbol(circle_hollow) || ///
scatter pptMean adj if Treated == 1, mcolor(navy) msymbol(circle)  ///
graphr(color(white)) legend( order( 5  "Untreated" 6 "Treated" 2 "Linear Fit, Unreated" 4 "Linear Fit, Treated" ) cols(2))  ///
xtitle("Adjudication Start Date") ///
ytitle("Mean Precipitation (mm)")



*Panel H: Temperature
twoway lfitci  tMean adj if Treated ==0, lcolor(maroon) lpattern(dash) alwidth(none) fcolor(gs13)       fintensity(inten50) || ///
lfitci tMean adj if Treated ==1, lcolor(maroon) alwidth(none) fcolor(gs13) || ///
scatter tMean adj if Treated == 0, mcolor(navy) msymbol(circle_hollow) || ///
scatter tMean adj if Treated == 1, mcolor(navy) msymbol(circle)  ///
graphr(color(white)) legend( order( 5  "Untreated" 6 "Treated" 2 "Linear Fit, Unreated" 4 "Linear Fit, Treated" ) cols(2))  ///
xtitle("Adjudication Start Date") ///
ytitle("Mean Temperature (deg C)")


restore


